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ABSTRACT 


The Planar Solid Oxide Fuel Cell (PSOFC) is one of the renewable energy 
technologies that is important as the main source for distributed generation 
and can play a significant role in the conventional electrical power 
generation. PSOFC stack modeling is performed in order to provide a 
platform for the optimal design of fuel cell systems. It is explained by the 
structure and operating principle of the PSOFC for the modeling purposes. 
PSOFC model can be developed using Artificial Neural Network approach. 


Keywords: The data required to train the neural net-work model is generated by 
ANN simulating the existing PSOFC model in the MATLAB/ Simulink software. 
The Radial Basis Function (RBF) and Multilayer Perceptron (MLP) neural 
Fuel cell networks are the most useful techniques in many applications and will be 
MLP applied in developing the PSOFC model. A detailed analysis is presented on 
PSOFC the best ANN network that gives the greatest results on the performances of 
RBF the PSOFC. The simulation results show that Multilayer Perceptron (MLP) 
gives the best outcomes of the PSOFC performance based on the smallest 

errors and good regression analysis. 
Copyright © 2019 Institute of Advanced Engineering and Science. 
All rights reserved. 
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1, INTRODUCTION 

The development of fuel cell as a renewable energy technology has become very popular due to its 
clean and sustainable sources of electrical energy and produces pure energy with comparatively high 
efficiency and effectively [1-3]. A fuel cell is a device that converts the chemical energy from a fuel such as 
methane into electrical energy for power applications through a chemical reaction of oxygen or any oxidizing 
agent with a positively charged hydrogen ions inside the cell [4]. It can be categorized as a first-rate 
electricity source because of its ability in producing continuous and constant power at full load. A fuel cell 
generation system needs oxygen and hydrogen in performing the chemical reaction inside the stack and 
produce electricity in DC form [5]. There are at least four fuel cell types have evolved in the past decades and 
becoming popular due to the clean and silent system. Fuel cell types are differentiate after their electrolyte, 
the substance that transports the ions and also dictates the operating temperature of a fuel cell types [6]. 
Depending on the operating temperature, a specific catalyst is chosen to oxidize the fuel such as Molten 
Carbonate Fuel Cell (MCFC), Phosphoric Acid Fuel Cell (PAFC), Solid Oxide Fuel Cell (SOFC) and Proton 
Exchange Membrane Fuel Cell (PEMFC). MCFC and SOFC operates at higher temperature and have the 
ability of internal fuel reformation that allows numerous fuel options. SOFC model can be divided into two 
main types, namely, tubular and planar configuration types. Tubular SOFC are formed into stacks by 
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implementing a tube sheet to support tubes, whereas the planar configuration is formed by a thin plate cells 
sandwiched between metal interconnecting plates [7]. Advantages of both types normally depend on their 
physical configurations. The tubular configuration is a matured technology which promises robustness, 
resistance to mechanical damage, and also has a simpler sealing, whereas the planar design is more compact 
and has a high power density, which decreases material content, thus lowering cost. Some tubular designs 
eliminate the need for seals and allow for thermal expansion. The tubular SOFC operates between 900°C- 
1000°C but planar operates at 800°C or even below and this main criteria effect the development cost of 
tubular SOFC which is normally higher [8, 9] and the thermal safety for tubular SOFC is quite challenging 
[10, 11]. In order to achieve high efficiency, planar SOFC normally consists of an ultra-thin electrode or 
electrolyte sheets with low electrical resistance. Operation at temperature less or lower than the common 
tubular SOFC enables less exotic materials of structure, thus cost saving [7]. 

There are many research works on developing tubular SOFC models such as a development of a 
dynamic model of SOFC for power system simulation using several types of software and also model for 
analyzing power system performances which consider all the losses or voltages drop in the model [12, 13]. 
Besides that, a nonlinear state space model of tubular SOFC to examine the dynamic behavior of SOFC for 
changes in input fuel, air pressure, and temperature was developed by Yutong et al. [14] and most of these 
researchers used MATLAB/Simulink environment in developing SOFC model due to its ability in analyzing 
the model for future power applications. The tubular types 1s receiving too much attention from researchers 
[15-17] while the modeling of planar SOFC is limited. 

PSOFC is the new focus in renewable energy technology where the PSOFC system balances the 
manufacturing cost and also promises the prospective in achieving lower stack. The main purpose of this 
paper is to model PSOFC performances by using artificial neural network (ANN) approach. Normally, in the 
existing technology, in order to develop the model of fuel cell, it will take a longer period of time to study the 
dynamic behavior due its complex system and the chemical reactions inside the model. Therefore, ANN can 
help in reducing the time required during the fuel cell modeling process. Moreover, it 1s basically difficult to 
estimate the parameters for the fuel cell system. Therefore, neural network approach will provide well 
modelling advantages, does not required many investigation of parameters and it will be easier in extending 
the model for power system analysis study such as distributed generation and micro grid applications. 


2. ELECTROCHEMICAL EQUATIONS FOR PLANAR SOFC 

Electrochemical equations are important in analyzing the PSOFC reactions inside the cell and all of 
these equations are used to estimate current densities throughout the cell. PSOFC Simulink model can be 
developed by using the electrochemical equations and the input output data for ANN model can be obtained 
from the simulation of the PSOFC Simulink model. PSOFC voltage is calculated based on the reversible 
voltage and the sum of several losses from the polarization curve [18] which are the effects of activation, 
con-centration and ohmic losses [19, 20]. The open circuit voltage (Eocy) or also known as the reversible cell 
voltage is acknowledged as a function of temperature and reactant concentration. The cell operating voltage 
(Vcei) 18 equal to the subtraction of the open circuit voltage from the sum of losses shown as below [18]: 


Veet = Eocv - >, (act + Yohm + Yconc) (1) 
The open circuit voltage can be obtained from the Nernst equation below [21]: 

Eocy=E° —RT/zF (In(Px2*Po2°°)/Px20) (2) 

where; E?= -AG?/zF (3) 


where AG° is the Gibbs free energy change at standard state and F refer to the Faraday’s number. n is the 
number of electrons formation in the reaction where the value of n 1s two due to the number of two electrons 
transferred for every one mole of reacted hydrogen. Normally in fuel cell characteristics as well as in SOFC 
studies, the major losses are the ohmic polarization (4onm), activation polarization (Yacr), and concentration 
polarization (Hconc). The activation loss is the difference between ionic potential at the electrolyte and 
electronic potential at the electrode which is greater than the equilibrium potential [22]. Activation 
polarization for both electrodes, anode and cathode can be determined using equation (4) and (5) below [22]: 


Nact ,an= 2RT/ZF [sinh! (I/2ioan )| (4) 


Nact. ca = 2RT/zF (sinh! (I2io.ca)] (5) 
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In the above equations, io.an and io,ca are the exchange current density for both anode and cathode. 


The ohmic polarization depends on the electrical resistance of the electrodes and the ion transport in the 
electrolyte and it is describe by the Ohm’s Law as shown below. 


oho = Roamb (6) 
The ohmic resistance, Rohm can be written as [22]: 


Ronn = [ (Tan/Can) +(Tei/Cel )+(Teca/Cca )] (7) 


Concentration polarization happens due to the reductions of reactant gas concentrations and can be 
determined by the following eequations [22]: 


Nconc,an=-(RT/2F )in[{ 1-U/ias)| +(RT/2 PF )ln[ 1+ (Px, an)/(Px20,an'tas) |] (8) 
Nconc,ca = -RT/4F In 1-W/ics)) (9) 


where igs and i;s are the current densities at electrode or electrolyte interfaces in which the partial pressure of 
gas 1s nearly zero and it can be define as [22]: 


ias=(2FP 32, an)(D (RT tan ) (10) 
ics=(4FP02,ca(D  \/([(Pea-Po2,0a) ))/Peal[ RT tea ] (11) 


The output fuel cell voltage, V;- can be expressed as [22]: 


Vic =No (Eocv- Nact,an-Y act,ca~4 ohm-Y conc,an-Y conc, ca) ( 1 2) 


The overall output voltage depends on the number of cells, (No) used in the design. Equations (1)- 
(12) are used to develop PSOFC model by using MATLAB/Simulink software as shown in Figure | and the 
parameters used in this developed model are shown in Table 1. The overall model of the PSOFC consists of 
seven subsystem which are partial pressure, open circuit voltage, activation polarization in anode, activation 
polarization in cathode, concentration polarization for both anode and cathode and also the ohmic 
polarization. The partial pressure equations can be found in El-Sharkh et.al [23]. 





moh 


Figure 1. PSOFC Simulation Model in MATLAB/ Simulink [24] 
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Table 1. PSOFC Model Parameters 


Parameters Values 
io.an (Am?) 5300 
inca (Am) 2000 
Tan (mm) 700 
Tei (mm ) 10 
Tca (mm ) 50 
Can (S/m ) (9x107/T)e!0 
Cet (S/m ) (3.34x104/T Je !9300/T 
Oca (S/m ) (4x107/T)e 1200/7 
Pea (atm ) 1 
Pan (atm ) 1 
5 
D gu Ge 3.5x10 
6 
‘5 ff (ms) 7.3x10 
Prp,an (Yomol ) 97 
Px20,an (Yomol ) 3 
Po2.ca (Yomol ) 21 
AG? 237kJmol! 


3. ANNINPUT AND OUTPUTS DATA 

An artificial neural network is a mapping network or technique that were able to compute some 
functional relationship between its input and output. Nowadays, the neural network have become a very 
importance tool or technique in some of the application as it were used to solve a wide variety of problems. 
Moreover, ANNs able to give results for non-linear problems very quickly once the training process is over 
[25-27]. Therefore, the performances of the PSOFC can be model by implementing the ANN approach 
employing MLP and RBF neural network techniques. The applications of the ANN is a relatively simple 
technique with a rapid response time and the system is considered suitable to determine the performances of 
the PSOFC. 

In the first part of the simulation, the PSOFC system have been design and model in the 
MATLAB/Simulink software in order to obtain the set of inputs and outputs for about 1200 data. These data 
later will be used for the next task which 1s to construct the PSOFC model using artificial neural network 
environment. The set of input and output data were divided into 3 parts which are the data of training, testing 
and validation. For ANN model, about 70% samples of data were used for training, 15% for testing and 15% 
for validation. For this paper, there are two inputs data involved which are the load current and temperature. 
The set of input data will be fed to the MLP neural network using sigmoid activation function in the hidden 
layer. There are two output neurons for the PSOFC system which are the output voltage and also the output 
power. Then number of neurons in hidden layer effects the overall neural network performances as small 
amount of neuron may have prohibited the network from correctly and effectively mapping the input and 
output. Table 2 shows the parameters used in MLP and RBF neural network and Figure 2 shows the basic 
arrangement of the MLP and RBF architecture. 


Table 2. MLP and RBF Parameters 


Methods MLP RBEF 
Goal (MSE) Not applicable 0 
Inputs 2 2 
Outputs 2 2 
Hidden layer 1 1 
Activation function in hidden layer Sigmoid Gaussian 
Spread Not applicable 8 
Training 840 840 
Testing 180 180 
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current voltage 
temperature—} — power 
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hidden layer 


Figure 2. The three layers architecture of MLP and REF neural networks 


Generally, normalization process is a process that transforming all the data or variable of a certain 
parameter to a specific range or interval like 0 and 1 or between -1 to +1. When there is big different in 
range, normalization is required so that all the inputs are at a comparable range. The data normalization is 
one of the important step for any analysis of the input sample from many situation in putting the data into a 
close range. In the ANN system, the existence of the activation function that used in the hidden layer with 
certain range made the data normalization process is one of the initial step before proceed to the next step. 
For sigmoid function, the range lies between 0 to | so the data will be normalized within the range in order to 
ensure the regression result lies within the specific area. Based on this paper, the normalization of data into 
p.u. from the range 0 until | can be done by using the equation below, 


Zi=(Xi-X (min) )/(X(max)-X(min)) (13) 


where, 
Zi = the normalized data 
X = the actual data 
The performance of the proposed ANN technique is evaluated in terms of the maximum error (max) 
and root mean square (RMS) error (@;ms) which are given mathematically by the following equations: 


€max= Max | Ty- Og |, q = 1,2,3....,n (14) 
€rm= V(1/n ¥ 4 -4 (tq-Ou)” ) (15) 


where Ty=[ tgi,tq2 ,.-- 5tamax ] 18 the target vector, Og=[Oq1,Oq2,.--,;Ogmax |] 18 the output vector of the output 
layer, and n is the number of neurons in the output layer. This error function is generated when the desired 
output of a network is compared with the actual output of a system. 


4. RESULTS AND ANALYSIS 

The results of the comparison between MLP and RBF for modeling the PSOFC performances are 
presented. In this section, the performance of the PSOFC system was determine by comparing the results of 
the MLP and RBF neural network model using a set of inputs and outputs data obtained from the PSOFC 
model simulation results using MATLAB/Simulink explained before. There are two main results that have 
been discussed which are the regression analysis and the network errors. The regression analysis is about the 
statistical process for estimating the relationship among the variables. This statistic will shows the 
distribution of data lies to the certain point in the regression process. The regression analysis of the MLP and 
RBF neural network between the outputs and responding targets are shown in Figures 3 and 4, respectively. 
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: R=0.99929 


Output ~= 1'Target + 0.00059 





Output ~= 0,96" Target + 0.027 





Figure 4. Regression of PSOFC model using RBF 


The correlation coefficient (R) for both ANN types are almost equal to 1 that indicates the 
regression have the strong correlation between the actual output data and target data of neural output. The 
MLP neural network indicates the best value of R which is 0.99929 followed by RBF function which is 
0.97893. From the regression results, it can observed that the MLP network employing sigmoid activation 
function shows the best result as the data are converges along the best fit line. The maximum and RMS errors 
are then calculated using equations (14) and (15) and Table 3 summarizes the network errors with bold 
figures represent the best testing RMS errors. 


Table 3. Performance of MLP and RBF Neural Networks Testing Stage 


OupubDat Sigmoid (MLP) Gaussian (RBF) 
Emax e rms Emax é rms 
Power 0.0692 0.0328 0.2822 0.0823 
Voltage 0.0623 0.0145 0.1021 0.0755 


From Table 3, it can be observed that the MLP network always achieves the best classification of the 
RMS errors for the two of the outputs data of PSOFC system. The MLP network is good in determining the 
testing stage for the performance of PSOFC system with the lowest RMS error. Based on the results, the 
MLP neural network shows the lowest value of the maximum error for power and voltage with the values of 


Indonesian J Elec Eng & Comp Sci, Vol. 15, No. 3, September 2019: 1645 - 1652 


Indonesian J Elec Eng & Comp Sci ISSN: 2502-4752 O 1651 


0.0692 p.u and 0.0623 p.u, respectively while for RBF neural networks, the maximum error for both power 
and voltage are 0.2822 p.u and 0.1021 p.u, respectively. MLP 1s also good in minimizing the maximum error 
for the testing stage with the value of 0.0328 p.u and 0.0623 p.u for both power and voltage compared to 
RBF with the higher values. This indicates that the MLP network with Sigmoid activation function in the 
hidden layer has the best results and able to give the best classification for the power and voltage values. It 
can be conclude that the trained MLP neural network can correctly predict the performances of PSOFC with 
the minimal errors. 


5. CONCLUSION 

The investigation of PSOFC performances by using ANN approach shows the performance of this 
fuel cell by comparing the results from different types of ANN which are the radial basis function and 
multilayer perceptron. Through the regression analysis it can be clearly seen the correlation amongst 
variables has been successfully estimated statically. It has been known practically that MLP is more superior 
than RBF but this is true for certain application as generically RBF is also powerful to yields best regression 
on given set data input and output depending how the training been conducted. To have minimum error for 
any network is utterly prominent thus it can be concluded through this investigation that MLP network does 
its best performance with showing less errors. In a nutshell, the MLP neural network with the Sigmoid 
activation function in the hidden layer gives the best result compared to the RBF as the regression shows the 
best fitting result with minimal RMS errors of 0.0692 p.u for power and 0.0623 p.u for voltage. 
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